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ABSTRACT. We establish a link between the study of completely integrable systems of partial 
differential equations and the study of generic submanifolds in C n . Using the recent developments 
of Cauchy-Riemann geometry we provide the set of symmetries of such a system with a Lie group 
structure. Finally we determine the precise upper bound of the dimension of this Lie group for 
some specific systems of partial differential equations. 
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1. Introduction 

To study the geometry of a real analytic Levi nondegenerate hypersurface M in C 2 , one of 
the principal ideas of H. Poincare, of B. Segre and of E. Cartan in the fundamental memoirs 
1201 . 1211 . 1221 . was to associate to M a system (£m) of (partial) differential equations, in 
order to solve the so-called equivalence problem. Establishing a natural correspondence between 
the local holomorphic automorphisms of M and the Lie symmetries of (£m) they could use the 
classification results on differential equations achieved by S. Lie in [5 1 and pursued by A. Tresse 
in l28l . 

Starting with such a correspondence, we shall establish a general link between the study of 
a real analytic generic submanifold of codimension m in C n+m and the study of completely 
integrable systems of analytic partial differential equations. We shall observe that the recent the- 
ories in Cauchy-Riemann (CR) geometry may be transposed to the setting of partial differential 
equations, providing some new information on their Lie symmetries. 

Indeed consider for K = M or C a K-analytic system (£) of the following general form: 

{£) u 3 xa (x) = F£ (x, u{x), {u 3 ^} q) (x))i<g< p ^ . 

Here x — {x\, . . . , x n ) £ IK™, u — (u 1 , . . . , u m ) £ K m , the integers j(l), . . . ,j(p) satisfy 
1 < j{q) < m for Q — L ■ ■ ■ ,P> an d a an d me multiindices (3(1), . . . ,(3(p) G N" satisfy 
M, \P(q)\ > 1. We also require (j, a) ^ (j(l), /3(1)), . . . , (j(p), For j = 1, . . . , m and 

a G N™, we denote by u 3 xa the partial derivative d^u 3 /dx a . We assume that the system (£ ) 
is completely integrable, namely that the Pfaffian system naturally associated in the jet space is 
involutive in the sense of Frobenius. We note that in that case (£) is locally solvable, meaning that 
through every point (x* , u* , u*p, u* a ) in thejet space, satisfying^* = F a (x* , u* , u*^) (writtenina 
condensed form), there exists a local K-analytic solution u = u(x) of {£) satisfying u(x*) = u* 
wdu x i}{x*) = u*p. Consequently the Lie theory (| 18 1) may be applied to such systems. We shall 
associate with (£) the submanifold of solutions M. in K™ +2m+p given by K-analytic equations 
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of the form 

(1) v? = Qj(x,u,x), j = l,...,m, 

where v 6 K OT and where x 6 K p . Moreover the integer m + p is the number of initial 
conditions for the general solution u(x) := tt(x,v,x) of (£), whose existence and unique- 
ness follow from complete integrability. Precisely, the parameters v : \ correspond to the data 
i<9<p- ^ n me special case where the system (£) is constructed from a generic 
submanifold M as in \2l\, |24| (see also Subsection 2.2 below), the corresponding submanifold 
of solutions is exactly the extrinsic complexification of M. 

A pointwise K-analytic transformation (x',vf) = $>(x,u) defined in a neighbourhood of 
the origin and sufficiently close to the iedntity mapping is called a Lie symmetry of (£) if it 
transforms the graph of every solution to the graph of an other local solution. A vector field 
X = Y^i=i Q l ( x ' u ) 9/dxi + Y^JLi RH X : u ) d/du^ is called an infinitesimal symmetry of {£) 
if for every s close to zero in K the local diffeomorphism (x, u) i— > exp(sX)(x, u) associated to 
the flow of X is a Lie symmetry of £ . According to [ 18 1 (Chapter 2) the infinitesimal symmetries 
of (£) form a Lie algebra of vector fields defined in a neighbourhood of the origin in K n x K m , 
denoted by 6rjm(£ ). Inspired by recent developments in CR geometry we shall provide in 
Section 2 nondegeneracy conditions on M insuring firstly that 6t)m(£ ) may be identified with 
the Lie algebra &ym(A4) of vector fields of the form 




which are tangent to Ai, and secondly that 6rjm(.M) = &t)m(£) is finite dimensional. The 
strength of this identification is to provide some (non optimal) bound on the dimension of 
&X)m(£ ) for arbitrary systems of partial differential equations with an arbitrary number of vari- 
ables, see Theorem[3] 

In the second part of the paper (Sections 3, 4 and 5), using the classical Lie theory (cf. l5l . 
1181 . [l9| and |2|), we provide an optimal upper bound on the dimension of 6rjm(£) for a 
completely integrable K-analytic system (£) of the following form: 

{£) U 3 x c = Fl(x,u(x), (Ua:/s(2O)l<|0|<re-l), « € N n , |«| = K, j = 1, . . . , TO. 

This system is a special case of the system studied in Section 2. For instance the homogeneous 
system (£q) : Ux k —xt, K ( x ) — is completely integrable. The solutions of (£ ) are the 
polynomials of the form u>(x) = X^gn™ \0\<k-i ^8 x>3 ' 3 ~ 1, . . . , m, where G K and a 
Lie symmetry of (£ ) is a transformation stabilizing the graphs of polynomials of degree < k — 1. 
We prove the following Theorem: 

Theorem 1. Let {£) be the WL-analytic system of partial differential equations of order k > 2, 
with n independent variables and m dependent variables, defined just above. Assume that (£) is 
completely integrable. Then the Lie algebra 6rjm(£ ) of its infinitesimal symmetries satisfies the 
following estimates: 

f dim K (6rim(£)) < (n + m + 2)(n + m), if K = 2, 

( ) \ dim K (6t)m(£)) < n 2 + 2i! + m 2 +mC„^_ 1 , if k > 3, 

where we denote C^] r l l _ 1 :— ^fj^zjjf- Moreover the inequalities become equalities for the 
homogeneous system (£q). 

We remark that there is no combinatorial formula interpolating these two estimates. Theorem^ 
is a generalization of the following results. For n = to = 1, S. Lie proved that the dimension 
of the Lie algebra 6rjm (£) is less than or equal to 8 if n = 2 and is less than or equal to n + 4 
if k > 3, these bounds being reached for the homogeneous system (cf. [5]). For n = 1, to > 1 
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and k = 2, F. Gonzalez-Gascon and A. Gonzalez-Lopez proved in 1 1 1 1 that the dimension of 
Srjm (£) is less than or equal to (m + 3)(m + 1). For n = 1, m > 1 and k = 2, using the 
equivalence method due to E. Cartan, M. Fels [6 1 proved that the dimension of 6rjm (£) is less 
than or equal to m 2 + 4m + 3, with equality if and only if the system {£) is equivalent to the 
system u° 2 = 0, j = 1, . . . , m. He also generalized this result to the case n = 1, m > 1, K = 3. 
For n. > I, m > 1 and k = 2, A. Sukhov proved in |24| that the dimension of 6rjm (£) is less 
than or equal to (n + m + 2)(n + m) (the first inequality in Theorem[0, with equality for the 
homogeneous system u J Xk Xk ^ = 0. 

Consequently, for the case k = 2, we will only give the general form of the Lie symmetries of 
the homogeneous system (£q) (see Subsection 5.2). We will prove Theorem[2for the case k > 3. 
The formulas obtained in Sections 3, 4 and 5 were checked with the help of MAPLE release 6. 

Acknowledgment. This article was written while the first author had a six months delegation 
position at the CNRS. He thanks this institution for providing him this research opportunity. The 
authors are indebted to Gerard Henry, the computer ingenieur (LATP, UMR 6632 CNRS), for his 
technical support. 

2. SUBMANIFOLD OF SOLUTIONS 

2.1. Preliminary. Let K = R or C. Let n > 1 and let x = (xi,. .., x n ) G N. We denote by 
K{x} the local ring of K-analytic functions ip = tp(x) defined in some neighbourhood of the 
origin in K n . If (p € K{x} we denote by ip the function in K{x} satisfying ip(x) = (p{x). Recall 
that a K-analytic function ip defined in a domain U C IK™ is called K-algebraic (in the sense 
of Nash) if there exists a nonzero polynomial P — P(X\, . . . , X n , $) € . . . , X n , $] 
such that P(x, f(x)) = on U. All the considerations in this paper will be local: functions, 
submanifolds and mappings will always be defined in a small connected neighbourhood of some 
point (most often the origin) in IK". 

2.2. System of partial differential equations associated to a generic submanifold of C" +m . 
Let M be a real algebraic or analytic local submanifold of codimension m in C n+m , passing 
through the origin. We assume that M is generic, namely TqM + iTqM = ToC n+m . Classically 
(cf. |1|) there exists a choice of complex linear coordinates t — (z,w) £ C™ x C m centered 
at the origin such that TqM = {Imaj ~ 0} and such that there exist m complex algebraic or 
analytic defining equations representing M as the set of (z, w) in a neighbourhood of the origin 
in C n+m which satisfy 

(4) wi = &i(z,z,w), ,w m = Q m (z,z,w). 

Furthermore, the mapping = (0i, . . . , ro ) satisfies the functional equation 

(5) w = 0(z, z, 6(z, z, w)), 

which reflects the reality of the generic submanifold M. It follows in particular from © that the 
local holomorphic mapping C m 9 w h- > (Oj(0, 0, w))i<j< m € C m is of rank m at w — 0. 

Generalizing an idea due to B. Segre in |21 1 and |22|, exploited by E. Cartan in 1 3 1 and more 
recently by A. Sukhov in 1 24 1 , 1 25 1 , 1 26 1 , we shall associate to M a system of partial differential 
equations. For this, we need some general nondegeneracy condition, which generalizes Levi 
nondegeneracy. Let 4fN with £q > 1. We shall assume that M is ^-finitely nondegenerate 
at the origin, cf. [\\, 1171 . 151 . This means that there exist multiindices . . . , (3(n) £ N n 
with \/3{k)\ > 1 for k — 1, . . . , n and maxi<fc<„ \(3{k)\ = £q, and integers j(l), . . . , j(n) with 
1 < j(k) < m for k = 1, . . . , n such that the local holomorphic mapping 

(6) c n+m 3 (z.w) ((e^o, z, w))i<j< m , (e jW .*««( ' 2 ' e £m+n 
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is of rank equal to n + m at (z,w) = (0,0). Here, we denote the partial derivative 
d'P'@j(0, z, w)/dz@ simply by 0^ z $ (0, z , w). Then M is Levi nondegenerate at the origin 
if and only if £q = 1. By complexifying the variables z and w, we get new independent vari- 
ables ( G C™ and £ G C m together with a complex algebraic or analytic m-codimensional 
submanifold M. in C 2 (™+ m ) of equations 

(7) Wj = Qj(z,C,£), j = l,...,m, 

called the extrinsic complexification of M. In the defining equations of AL following [21 1 
and 1 24 1, we may consider the "dependent variables" w\, . . . , w m as algebraic or analytic func- 
tions of the "independent variables" z = (z\, ■ ■ ■ , z n ), with additional dependence on the ex- 
tra "parameters" (C, £) G <C n+m . Then by applying the differential operator d^/dz a to Q, 
we obtain Wj >z a (z) = Qj !Z a(z,£,Q- Writing these equations for (j, a) = (j(k), (3(k)) with 
k = 1, . . . , n, we obtain a system of rn + n equations 

f Wj(z) = &j(z,C,€), j = l,...,m, 

\ Wj{k),zPV*){z) = ®j{k)^*){ z X,t), k = l,...,n. 

In this system (|8j, by the assumption of ^o-finite nondegeneracy 0, the algebraic or ana- 
lytic implicit function theorem allows to solve the parameters (£, £) in terms of the variables 
(zk, Wj(z), Wjt k -\ z p{k) (z)), providing a local algebraic or analytic C™ +m -valued mapping R 
such that (C;0 = R( z k,Wj(z),Wj^ k ^ z i3(k)(z)y Finally, for every pair (j,a) different from 
(1,0), . . . , (m,0), (j(l),/3(l)), • ■ ■ , {j(n),l3{n)) 7 we may replace (CO by R in the differenti- 
ated expression Wj. z a (z) — Qj. z ^ (z, ^). This yields 

w Ji2Q (z) = 9j,^ (z, J R(z fc ,u> i (z),w i(fe)z{i ( fc )(z))) 

=: Fj,a (zk,Wj(z),Wj( k ), z 0V>)(z)) ■ 

This is the system of partial differential equations associated with Ai. As argued by B. Segre 
in 1211 . the geometric study of generic submanifolds of C n may gain much information from 
the study of their associated systems of partial differential equations (cf. |24|, |25|). The next 
paragraphs are devoted to provide a general one-to-one correspondence between completely 
integrable systems of analytic partial differential equations and their associated "submanifolds 
of solutions" (to be defined precisely below) like M. above. Afterwards, we shall observe that 
conversely, the study of systems of analytic partial differential equations also gains much infor- 
mation from the direct study of their associated submanifolds of solutions. 

2.3. Completely integrable systems of partial differential equations. Let now n, m, p G N 

with n, m, p > 1, let n G N with n > 2 and let u — (u 1 , . . . , u m ) G K m . Consider a collection 
ofpmultiindices /3(1), . . . ,(3(p) G N™ with \/3(q)\ > lforg = 1, . . . ,pandmaxi< g < p \/3(q)\ = 
K — 1. Consider also p integers . . . , j(p) with 1 < j(q) < m for q — 1, . . . ,p. Inspired 
by (|9jl, we consider a general system of partial differential equations of n independent variables 
(xi, . . . , x n ) and m dependent variables (u 1 , . . . , u m ) which is of the following form: 

(£) U J x a (x) = Fl (x, U(x), (u?J*h (x))l<q< p ) , 

where (j, a) ^ (3(1)), . . . , (j(p), (3(p)) and j = 1, ...,m, \a\ < K. Here, we assume 

that u = is a local solution of the system (£) and that the functions are K-algebraic or 
IK-analytic in a neighbourhood of the origin in K n + m +P. Among such systems are included 
ordinary differential equations of any order k > 2, systems of second order partial differential 
equation as studied in |24|, etc. 

Throughout this article, we shall assume the system (£) completely integrable. By analyzing 
the application of the Frobenius theorem in jet spaces, one can show (we will not develop this) 
that the general solution of the system (£ ) is given by u(x) :— fl(x, v, x), where the parameters 
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v G K™ andx € K n essentially correspond to the "initial conditions" u(0) and (u 3 '^' a3 (0))k 9 < p , 
and ft is a K-analytic K7 l -valued mapping. In the case of a generic submanifold as in Subsec- 
tion 2.2 above, we recover the mapping O. In the sequel, we shall use the following terminology: 
the coordinates (x, u) will be called the variables and the coordinates (is, x) will t> e called the 
parameters or the initial conditions. In Subsection 2.5 below, we shall introduce a certain duality 
where the roles between variables and parameters are exchanged. 

2.4. Associated submanifold of solutions. The existence of the function fl and the analogy 
with Subsection 2.2 leads us to introduce the submanifold of solutions associated to the com- 
pletely integrable system (£), which by definition is the m-codimensional K-analytic submani- 
fold of K" +2m+p , equipped with the coordinates (x, u, is, x), defined by the Cartesian equations 

(10) Uj =n j (x,v,x) ) j = l,...,m. 

Let us denote this submanifold by A4. We stress that in general such a submanifold cannot 
coincide with the complexification of a generic submanifold of C" l+n , for instance because K 
may be equal to K or, if K = C, because the integer p is not necessarily equal to n. Also, 
even if IK = C and n = p, the mapping £1 does not satisfy a functional equation like In 
fact, it may be easily established that the submanifold of solutions of a completely integrable 
system of partial differential equations like (£) coincides with the complexification of a generic 
submanifold if and only if K = C, p = n and the mapping il satisfies a functional equation 
like 0. 

Let now M. be a submanifold of K n + 2 "+P G f the form i ll Oi l, but not necessarily constructed as 
the submanifold of solutions of a system (£). We shall always assume that flj (0, is, x) = v 3 ■ We 
say that Ai is solvable with respect to the parameters if there exist multiindices (3( 1 ),..., (3(p) <E 
N ra with \(3(q)\ > 1 for q = 1, . . . ,p and integers j(l), . . . , j(p) with 1 < j(q) < m for 
q = 1, . . . ,p such that the local K-analytic mapping 

(11) K m+ P 3 (is, x) — ► ((% (0, v, X )i<i< m , (%,),*«-> (0, v, X )) G K m +r 

is of rank equal to m + p at (£, x) — (0, 0) (notice that since f2j(0, v, x) = vK then the m first 
components of the mapping (lilt are already of rank m). We remark that the submanifold of 
solutions of a system (£) is automatically solvable with respect to the variables, the multiindices 
(3(q) and the integers j(q) being the same as in the arguments of the right hand side terms in 
{£). 

2.5. Dual system of defining equations. Since £lj (0, v, %) = , we may solve the equa- 
tions dlOl with respect to v by means of the analytic implicit function theorem, getting an equiv- 
alent system of equations for A4: 

(12) v 3 ' = tt*j(x,x,u), j = l,...,m. 

We call this the dual system of defining equations for M. . By construction, we have the functional 
equation 

(13) u = n(x,n*(x,x,u),x), 

implying the identity SI* (0, x, u) = uK We say that M. is solvable with respect to the variables 
if there exist multiindices 6(1), . . . , 5(n) G N p with \8(l)\ > 1 for / = 1, . . . , n and integers 
j(l), ■ ■ ■ , j(n) with 1 < j(l) < m for I = 1, . . . , m such that the local K-analytic mapping 

(14) K"+ m 3 (X,U) .— » ( (Cllffi, X^lKjKrn, x ,(0 (0, X, tt)) ^ J £ K m+n 

is of rank equal to n + m at (x, u) = (0, 0) (notice that since f2*(0, x, u) = u- 7 , the m fisrt 
components of the mapping dl4> are already of rank m). 
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In the case where M. is the complexification of a generic submanifold then the solvability with 
respect to the parameters is equivalent to the solvability with respect to the variables since fi* = 
fi. However we notice that a submanifold M of solutions of a system (£) is not automatically 
solvable with respect to the variables, as shows the following trivial example. 

Example 1. Let n = 2, to = 1 and let (£) denote the system u X2 — 0, u XlXl — 0, whose general 
solutions are u(x) = v + x\\ ='■ Vl{x\,xi, v, \). Notice that the variable X2 is absent from the 
dual equation v = u — xixi =■ ^*(X; x i> x 2, u). It follows that A4 is not solvable with respect 
to the variables. 

2.6. Symmetries of (£), their lift to the jet space and their lift to the parameter space. We 

denote by J£ m the space of jets of order n of K-analytic mappings u = u{x) from K" to K' Tl . 
Let 

(15) {x h u*,Uil,Uil h , J G JK"+ mC " + " 

denote the natural coordinates on J£ m - Here, the superscripts j,ii and the subscripts 
I, ?i, ?2, • • • ; Ik satisfy j, i\ — 1, . . . , m and l,h,l2, . • ■ ,l K — 1, . . . , n. The independent co- 
ordinate lx corresponds to the partial derivative Xl . Finally, by symmetry of partial 

differentiation, we identity every coordinate ^ with the coordinates ^(Zi) <?{i\Y wnere 
a is an arbitrary permutation of the set {1, . . . , A}. With these identifications, the K-th order jet 
space Jn,m i s °f dimension n + m C£ +n , where Cj> :— q , (^'_ q y_ denotes the binomial coefficient. 
Also, we shall sometimes use an equivalent notation for coordinates on J£ m '- 

(16) ( IIl tf,[/^K" +mC .V i 

where (3 G N™ satisfies \(3\ < k and where the independent coordinate Up corresponds to the 
partial derivative v? . 

associated to the system (£) is the so-called skeleton Ag, which is the K-analytic submanifold 
of dimension n + m + p in J^ m simply defined by replacing the partial derivatives of the 
dependent variables v? by the independent jet variables in (£): 

(17) Ui=Fi(x,u,(U$ ) ) 1 < q < p ), 

for (j, a) f (j(l), . . . , (j(p), (3{p)) and j = 1, . . . , to, |a| < k. Cleai'ly, the natural 

coordinates on the submanifold of J^m are tne n + m + p coordinates 

(is) (x^,(u;{i\) 

Let h — h(x,u) be a local K-analytic diffeomorphism of K n+m close to the identity mapping 
and let ir K : J^ m — » K rl+m be the canonical projection. According to 1 18 1 (Chapter 2) there 
exists a unique lift of h to J£ m such that 7r K o —ho n K . The components of 
may be computed by means of universal combinatorial formulas and they are rational functions 
of the jet variables dl5> . their coefficients being partial derivatives of the components of h, see 
for instance §3.3.5 of 1 2 1 . By definition, h is a local symmetry of (£) if h transforms the graph of 
every local solution of (£) into the graph of another local solution of (£). This definition seems 
to be rather uneasy to handle, because of the abstract quantification of "every local solution", but 
we have the following concrete characterization for h to be a local symmetry of (£), cf. Chapter 2 
in fT8l. 

Lemma 1. The following conditions are equivalent: 

(1) The local transformation h is a local symmetry of (£). 

(2) Its K-th prolongation h^ is a local self-transformation of the skeleton Ag of{£). 
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These considerations have an infinitesimal version. Indeed, let X — Y2i=i Q l ( x i u ) d/dxi + 
2j=i R?{ x , u ) d/du^ be a local vector field with K-analytic coefficients which is defined in 
a neighbourhood of the origin in K n+m . Let s £ K and consider the flow of L as the one- 
parameter family h s (x, u) := exp(s X)(x, u) of local transformations. We recall that X is an 
infinitesimal symmetry of (£) if for every small s £ K, the mapping h a (x, u) := exp(s X)(x, u) 
is a local symmetry of {£). By differentiating with respect to s the K-th prolongation (h s )^ 
of h s at s = 0, we obtain a unique vector field X( K ) on the K-th jet space, called the n-th 
prolongation of X and which satisfies (-Kk)*{X^) = X. In Subsections 3.1 and 3.2 below, we 
shall analyze the combinatorial formulas for the coefficients of X^ K \ since they will be needed 
to prove Theorem^ 

Let Xg be the projection to the restricted jet space K m + n +?, equipped with the coordi- 
nates i ll 8t . of the restriction of X^ to As, namely 

(19) X £ := {n K , p )*(X^\ A£ ). 

The following Lemma, called the Lie criterion, is the concrete characterization for X to be an 
infinitesimal symmetry of (£) and is a direct corollary of Lemma^ cf. Chapter 2 in 1 18]. This 
criterion will be central in the next Sections 3, 4 and 5. 

Lemma 2. The following conditions are equivalent: 

(1) The vector field X is an infinitesimal symmetry of {£). 

(2) Its n-th prolongation X^ is tangent to the skeleton Ag. 

We denote by 6rjm(£) the set of infinitesimal symmetries of (£). Since it may be easily 
checked that (cX + dY)M = cX^ + dY^ and that [X^\ F (k) ] = ([X, F]) (k) , see Theo- 
rem 2.39 in 1 18 1, it follows from Lemma[2](2) that Grjm(£) is a Lie algebra of locally defined 
vector fields. Our main question in this section is the following: under which natural conditions 
is &X}m{£) finite- dimensional ? 

Example 2. We observe that the Lie algebra Srjm(£) of the system (£) presented in Example^ 
is infinite-dimensional, since it includes all vector fields of the form X = Q 2 (xi, x-x, u) d/dx 2 , 
as may be verified. As we will argue in Proposition [2] below, this phenomenon is typical, the 
main reason lying in the first order relation u X2 = 0. 

By analyzing the construction of the submanifold of solutions M. associated to the system 
(£), we may establish the following correspondence (we shall not develop its proof). 

Proposition 1. To every infinitesimal symmetry X = $^I==l Q l ( x ,u) d/dxi + 
Y)j—i R^(x, u) d/du^ of {£), there corresponds a unique vector field of the form 

m f) p f) 

(20) X = E^^ X ) — + J2A^, X )—, 

3=1 9=1 Aq 

whose coefficients depend only on the parameters {y, x)> such that X + X is tangent to the 
submanifold of solutions A4. 

This leads us to define the Lie algebra Srjm(TW) of vector fields of the form 

(2d E^m^+E Rj «) ^ + £ ff x)^-+E A x) w 

1=1 j=l 3 = 1 q=l 

which are tangent to M.. We shall say that the submanifold Ai is degenerate if there ex- 
ists a nonzero vector field of the form X = Ya=i Q l {x,u)d/dxi + Y^JLi R?( x ,u) d/dv? 
which is tangent to AL which means that the corresponding X part is zero. In this case, we 
claim that Srjm(.M) is infinite dimensional. Indeed there exists then a nonzero vector field 
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T = Y^iQ ( x > u )9/dxi + Y^ijLi IU{ x i u)d/du 3 tangent to M.. Consequently, for every 
K-analytic function A(x, u), the vector field A(x, u) T belongs to &t]m(A4), hence ©rjm(A^) 
is infinite dimensional. 

By developing the dual defining functions of M. with respect to the powers of \, we may 
write 

(22) = Sl*(x, x,u)=J2 X 7 ^* 7 (^ «), 

where the functions f2* (x, u) are K-analytic in a neighbourhood of the origin, we may formu- 
late a criterion for M. to be non degenerate with respect to the variables (whose proof is skipped). 

Proposition 2. The submanifold M is not degenerate with respect to the variables if and only if 
there exists an integer k such that the generic rank of the local ^-analytic mapping 

W (*>") » ( n h( X > U ))l<j< m , 7 GNP, | 7 |<fc 

is equal to n + m. 

Seeking for conditions which insure that &t)m(Ai) is finite-dimensional, it is therefore natural 
to assume that the generic rank of the mapping (1231 is equal to n + m. Furthermore, to simplify 
the presentation, we shall assume that the rank at (x, u) = (0, 0) (not only the generic rank) of 
the mapping H231 is equal to n + m for k large enough. This is a "Zariski-generic" assumption. 
Coming back to 14l >. we observe that this means exactly that M is solvable with respect to the 
variables. Then we denote by £q the smallest integer k such that the rank at (x, u) = (0, 0) 
of the mapping (I23> is equal to n + m and we say that M. is £y-solvable with respect to the 
variables. Also, we denote by £q the integer maxi< 9 < p \0(q)\ and we say that M. is £o-solvable 
with respect to the parameters. 

2.7. Fundamental isomorphism between Srjm(£) and ©lim(7W). In the remainder of this 
Section 2, we shall assume that M is £o-solvable with respect to the parameters and ^g-solvable 
with respect to the variabes. In this case, viewing the variables (y l , . . . , v m ) in the dual equations 
vi = ilj(x, x, u) of M. as a mapping of x with (dual) "parameters" (x, u) and proceeding as 
in Subsection 2.2, we may construct a dual system of completely integrable partial differential 
equations of the form 

(£*) 4-(x) = Gir (x. "(X), (< ( ,?o (x))i<Kn) , 

where (j, 7) ^ U0-), ^(1)), • • • , 3( n ))- This system has its own infinitesimal symmetry 

Lie algebra 6nm(£ *). 

Theorem 2. If M. is both solvable with respect to the parameters and solvable with respect to 
the variables, we have the following two isomorphisms: 

(24) 6tjm(£) = 6tjm(.M) = ©tjm(f*), 

namely X < > X + X < > X. 

In Subsection 2. 10 below, we shall introduce a second geometric condition which is in general 
necessary for Gr)m(.M) to be finite-dimensional. 

2.8. Local (pseudo)group Sym(A / () of point transformations of M. We shall study the ge- 
ometry of a local K-analytic submanifold M. of ]K™+ 2m +P whose equations and dual equations 
are of the form 



(25) 
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Let t := (x, u) G K" +m and r := {v, x) G K" +m . We are interested in describing the set of 
local K-analytic transformations of the space JC™+ 2m +J> which are of the specific form 

(26) (*,t).— 

and which stabilize .M, in a neighborhood of the origin. We denote the local Lie pseudogroup 
of such transformations (possibly infinite-dimensional) by Sym(.A4). Importantly, each transfor- 
mation of Sym(.M) stabilize both the sets {t — ct.} and the sets {r = ct.}. Of course, the Lie 
algebra of Sym(M) coincides with 6rjm(.A4) defined above. 

2.9. Fundamental pair of foliations on M. Letp G K" +2m+p be a fixed point of coordinates 
(t po ,t po ). Firstly, we observe that the intersection M. n {r = r po } consists of the n-dimensional 
K-analytic submanifold of equation u = Q(x,t po ). As r po varies, we obtain a local K-analytic 
foliation of M by n-dimensional submanifolds. Let us denote this first foliation by T v and 
call it the foliation of M. with respect to parameters. Secondly, and dually, we observe that the 
intersection M fl {t — t Po } consists of the p-dimensional K-analytic submanifold of equation 
v = 0*(x, t Po ). As t Po varies, we obtain a local K-analytic foliation of M by p-dimensional 
submanifolds. Let us denote this second foliation by T v and call it the foliation of M. with 
respect to the variables. We call {Fp, T v ) the fundamental pair of foliations on M.. 

2.10. Covering property of the fundamental pair of foliations. We wish to formulate a geo- 
metric condition which says that starting from the origin in M and following alternately the 
leaves of T v and the leaves of T v , we cover a neighborhood of the origin in Ai . Let us introduce 
two collections {Ck)\<k<n and (C* q )x< q < p of vector fields whose integral manifolds coincide 
with the leaves of T p and T v : 

d dflj , . d 

(27) I 3 

a ™ dn* d 

Let po be a fixed point in M of coordinates (x Po , u PQ , v PQ , \va ) e ^>i+2m+ P ^ j et Xi ._ 

(xi.i, . . . , x\ t n) G K n be a "multitime" parameter and define the multiple flow map 

(28) 

£ Xl (x Pa ,u PQ ,is Po ,Xp ) ■= exp(xi£)(p ) :=exp(a;i 1 „£ n (---(exp(a;i,i£i(po)))---)) : = 
:= (x Po + xi, Q(x Po + xi, v Po , Xpo)i v vo > Xpo)- 
Similarly, for x = (xi,i> ■ • ■ > Xi,p) € IK P ' define the multiple flow map 

(29) , Up Q , fp , Xpo ) '~ ( X P0 J U P0 5 

We may define now the mappings which correspond to start from the origin and to move alter- 
nately along the two foliations T p and T v . If the first movement consists in moving along the 
foliation T v , we define 

ri(a;i) :=C X1 (0), 

T.ixuxi) :=C* Xl (C Xl (0)), 
r 3 (xi,xi,a;2) := C X2 {C* xl (C Xl (ti))), 
r 4 (x 1; xi, x 2 , X2) ■= c* X2 (C X2 (C X1 (C X1 (0)))). 

Generally, we may define the maps Tk([xx]k), where [xx]k — Xi> x 2, X2, ■ • ■ ) with exactly 
k terms and where each xi belongs to K n and each xi belongs to W. On the other hand, if 
the first movement consists in moving along the foliation T p , we start with r^(xi) := C Xi (0), 
r^xij^i) := jC Xi (£ Xi (0)), etc., and generally we may define the maps Tl([xx]k), where 



k = 1, . . . , n, 
k = 1, . . . , n. 
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[x x ]k = (xii x i> X2j X2, ■ ■ ■ ), with exactly k terms. The range of both maps Tk and T^ is 
contained in M. We call Tk the k-th chain and T£ the k-th dual chain. 

Definition 1. The pair of foliations (T P ,J- V ) is called covering at the origin if there exists an 
integer k such that the generic rank of Tk is (maximal possible) equal to diniK Ai. Since the 
dual (fc + l)-th chain for xi = identifies with the fc-th chain Tk, it follows that the same 
property holds for the dual chains. 

In terms of Sussmann's approach |27|, this means that the local orbit of the two systems of 
vector fields (£fc)i<fc<n and (£q)i<q<p is °f maximal dimension. Reasoning as in |27| (using 
the so-called backward trick in Control Theory, see also 1171 1. it may be shown that there exists 
the smallest even integer such that the ranks of the two maps T2 Po and rjjL at the origin 
(not only their generic rank) in K nA10+PAI ° are both equal to dimjK M.. This means that r2 Mo and 
T^ are submersive onto a neighborhood of the origin in M. We call [iq the type of the pair of 
foliations [T v , T v ). It may also be established that fio < m + 2. 

Example 2.46. We give an example of a submanifold which is both 1-solvable with respect to 
the parameters and with respect to the variables but whose pair of foliations is not covering: with 
n = 1, m = 2 and p = 1, this is given by the two equations u 1 = v 1 , u 2 = v 2 + x\i- Then 
&t)m(A4) is infinite-dimensional since it contains the vector fields a(u 1 ) d/du 1 + a(i/ 1 ) d/dv 1 , 
where a is an arbitrary IK-analytic function. For this reason, we shall assume in the sequel that 
the pair of foliations {T p , T v ) is covering at the origin. 

2.1 1. Estimate on the dimension of the local symmetry group of the submanifold of solu- 
tions. We may now formulate the main theorem of this section, which shows that, under suitable 
nondegeneracy conditions, Sym(A^) is a finite dimensional local Lie group of local transforma- 
tions. If t e K n+m , we denote by \t\ := max 1 <j t < n+m \t^\. If (h, (f>) £ Syra(M) we denote 
by J t fc /i(0) the fc-th order jet of h at the origin and by J^0(O) the fc-th order jet of <f) at the origin. 
Also, we shall assume that M. is either K-algebraic or IK-analytic. Of course, the IK-algebraicity 
of the submanifold of solutions does not follow from the K-algebraicity of the right hand sides 
Fl of the system of partial differential equations {£). 

Theorem 3. Assume that the K-algebraic or ^-analytic submanifold of solutions A4 of the 
completely integrable system of partial differential equations {£) is both l^-sovable with respect 
to the parameters and i^-solvable with respect to the variables. Assume that the fundamental 
pair of foliations (J^p,^) is covering at the origin and let fiQ be its type at the origin. Then 
there exists So > such that for every e with < e < £q, the following four properties hold: 

(a) The (pseudo)group Sym(.M) of local ^-analytic diffeomorphisms defined for {(t, r) G 
K n+2m+ P . | t | < £j | T | < £ j w j qicn are f the form (t,r) h-> (h(t) , 4>(t)) and which 
stabilize M. is a local Lie pseudogroup of transformations of finite dimension del 

(b) Let kq '■= fJ-oi^o + ^o)- Then there exist two ^L-algebraic or ^-analytic mappings H Ko 
and $ Ko which depend only on A4 and which may be constructed algorithmically by 
means of the defining equations of M. such that every element (h, <fi) € Sym(A^), suffi- 
ciently close to the identity mapping, may be represented by 



Consequently, every element of Sym(A^) is uniquely determined by its K -th jet at the 
origin and the dimension d of the Lie algebra 6rjm(A^) is bounded by the number of 
components of the vector ( J 4 K °/i(0), J?° (f)(0)), namely we have 



(31) 



( 



h(t)=H K0 (t, J t Ko h(0)), 
^(r)=$ K0 (r, J£V(0)). 



(32) 



dim K 6rjm(£) = dim K &t)m(M) < (n + m) C\ 



1KQ 



+ (to + p) C\ 



1KQ 
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(c) There exists e' with < e' < e and a W-algebraic or ¥L-analytic mapping (Hm , <i>_A4 ) 
which may be constructed algorithmically by means of the defining equations of M, 
defined in a neighbourhood of the origin in ]g n + 2m +P x K d with values in ]K n+2m+p 
and which satifies (i/y^i, 0), $x(r, 0)) = (t,r), such that every element (h,4>) G 
Sym(A^) defined on the set {(i,r) G K n+2m+ P : \t\ < e' , \t\ < s'}, sufficiently 
close to the identity mapping and stabilizing M. may be represented as (h(t), 0(t)) = 
(H_M(t, St <p), &m( t , Sh cf>)) far a unique element Sh <f> G K. d depending on the mapping 
(h,6). 

(d) The mapping (t,T,s) i — > (Hm (t, s), $x(r, s)) defines a local K-algebraic or K- 
analytic Lie group of local ^-algebraic or ^-analytic transformations stabilizing M.. 

2.12. Applications. The proof of Theorem[3] which possesses strong similarities with the proof 
of Theorem 4.1 in |8 1, will not be presented. It seems that Theorem[3] together with the argu- 
mentation on the necessity of assumptions that M be solvable with respect to the variables and 
that its fundamental pair of foliations be covering, is a new result about the finite-dimensionality 
of a completely integrable system of partial differential equations having an arbitrary number of 
independent and dependent variables. The main interest lies in the fact that we obtain the algo- 
rithmically constructible representation formula (13 1 1 together with the local Lie group structure 
mapping (Hm, ®m)- In particular, we get as a corollary that every transformation (h(t), 4>(t)) 
given by a formal power series (not necessarily convergent) is as smooth as the applications 
(iJ Ko , 3> Ko ) are, namely every formal element of Sym(A^) is necessarily K-algebraic or K- 
analytic. As a counterpart of its generality, Theorem 3 does not provide optimal bounds, as 
shows the following illustration. 

Example 2.46. Let n = m = 1, let K > 3 and let {£) denote the ordinary differential equation 
u xK (x) = F(x,u(x),u x (x), . . . ,u xK -i(x)). Then the submanifold of solutions Ai is of the 
form u — v + xxi + ■ ■ ■ + x K, ~ 1 x K -i + 0(|a;| K ) + 0(|x| 2 )- It may be checked that Iq = k — 1, 
£q = 1 and /io = 3, hence kq — 3k. Then the dimension estimate in d32i is: dimjj 6rjm(£ ) < 
2 C|" 3k + k Cf£. This bound is much larger than the optimal bound dimK Grjm(£) < K + 4 
due to S. Lie (cf. [5 1; see also the case n = m = 1 of Theorem^. 

Untill now we focused on providing the set of Lie symmetries of a general system of partial 
differential equations with a local Lie group structure. As a byproduct we obtained the (non 
optimal) dimensional upper bound ( 1321 of Theorem[3] In the next Sections 3, 4 and 5, using the 
classical Lie algorithm based on the Lie criterion (see Lemma|2j, we provide an optimal bound 
for some specific systems of partial differential equations, answering an open problem raised 
in GD page 206. 

3. Lie theory for partial differential equations 

3.1. Prolongation of vector fields to the jet spaces. Consider the following K-analytic system 
(£) of non linear partial differential equations: 

( 33 ) u x kl -x kK i x ) = F L,-,kn ( x ' M ( a: )'< 1 !l ( a: )'---'< 1 il -x Ire _ 1 ( iE )) . 

where 1 < ki < ■ ■ ■ < k K < n, 1 < j < m, and Ff. k are analytic functions of 

n + ™C r ^+K-i variables, defined in a neighbourhood of the origin. We assume that (£) 
is completely integrable. The Lie theory consists in studying the infinitesimal symmetries 
X = 2~^/li Q l (x,u)d/dxi + Y^jLi {x , u) d / dui of (£). Consider the skeleton of (£), 
namely the complex subvariety of codimension m C% + n _j in the jet space J£ m > defined by 



(34) 
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where = 1, . . . , m and k\, . . . , k K , h, . . . , l K -i = 1, . . . , n. For k = 1, . . . , n let D k be the 
fc-th operator of total differentiation, characterized by the property that for every integer A > 2 
and for every analytic function P = P(x, u, Uf 1 , . . . , UJ 1 t ) defined in the jet space ^T^rn' 
the operator D k is the unique formal infinite differential operator satisfying the relation 



(35) 



' [D k P] (*),•••,<• 



(x) 



d 

dx k 



P 



(x) 



Note that this identity involves only the troncature of Dk to order A, denoted by D^, and defined 
by 



(36) 



i± — 1 ?i— 1 /i — 1 



E E 

ii = l ! 1 ,...,i A _ 1 =l 



K,l\,...,l\— i QJjil 



d 



According to Theorem 2.36 of IT51 . the prolongation of order n of a vector field X = 
Y^i—i Q l ( x : u ) d/dxi + Y^ijLi u) d/du J , denoted by X( K \ is the unique vector field on 

the space 3£ m of the form 



(37) 



m 71 771 71 r, 

^'^EE^i+E E Km " 

ki j=l k!,k 2 = l 



3 = 1 k 1= l 



+ 



E E R ki,...,k K 

j = l fei fc K = l 



d 



^k!,k 2 ,...,k K 



corresponding to the infinitesimal action of the flow of X on the jets of order k of the graphs of 
maps u = u(x), and whose coefficients are computed recursively by the formulas 



(38) 



h=i 

71 

D i 2 ( R iJ - E D L(Q h ) u L,h> 



R il,fe: 



R fei,fc 3 ,...,^ '— D kA R k u -M-^ ~ E D lx(Q h ) U k,~,kx-i,h- 

h=l 

For a better comprehension of the general computation, let us start by computing R K in the case 
n = m = 1. 

3.2. Computation of R K when n = m = 1. A direct application of the preceding formulas 
leads to the following classical expressions: 

(R 1 =R X + [R u - Q x ] U 1 + \-Q u ] {Uxf. 

(39) J R 2 = R X 2 + [2R XU - Q X 2] U 1 + [R u 2 - 2Q XU ] (U 1 ) 2 + [-Q u *] (U 1 )'^ 

{ + [R U -2Q X ]U 2 + [-ZQ u ]U l U 2 . 

Observe that these expressions are polynomial in the jet variables, their coefficients being differ- 
ential expressions involving a partial derivative of R (with a positive integer coefficient) and a 



SYMMETRIES OF PARTIAL DIFFERENTIAL EQUATIONS 



13 



partial derivative of Q (with a negative integer coefficient). We have also: 

' R 3 = R x s + [3R X 2 U - QxA U 1 + [3R XU 2 - 3Q X 2 U ] (U v f+ 

+ [R u s - 3Q XU 2] (U 1 ) 3 + [-Q u8 ] (U 1 ) 4 + [3R XU - 3Q X *} U 2 + 
+ [3R U 2 - 9Q XU ] U 1 !/ 2 + [-6Q U 2] (U 1 ) 2 !/ 2 + [-3Q a ] (U 2 ) 2 + 
+ [R u - 3Q X ] U 3 + [-4QJ U X U 3 . 
R 4 = R X 4 + [AR x s u ~ Q x ,\ U 1 + [&R X 2 U 2 - 4Q x3u ] (U 1 ) 2 + 
(40) { +[4,R xu 3-6Q x 2 u 2}(U 1 ) 3 + [R u .-4Q xu 3](U 1 ) 4 + [-Q ui ](U 1 ) 5 + 

+ [6R X 2 U - 4Q x a] U 2 + [12R xu i - 18Q X 2 U ] U 1 U 2 + 
+ [6R u3 - 2AQ XU 2] (U^U 2 + [-10Q&] {U 1 ) 3 U 2 + 
+ [3R U 2 - \2Q XU ] (U 2 ) 2 + [-15Q„2] U\U 2 ) 2 + [AR XU - 6Q x2 ] U 3 + 
+ [AR U 2 - 16Q XU ] U 1 ^ 3 + [-1QQ V *] (U^U 3 + [-10Q„] U 2 U 3 + 
+ [Ru - 4QJ U 4 + [-5Q U ] V^U 4 . 

Remark that all the brackets involved in equations J40i are of the form [A R x a u b+i — fi Q x a+i u b], 
where A,fi6N and a, b e N. 

In what follows we will not need the complete form of R K but only the following partial form: 

Lemma 3. For k > 4: 

K K = R x k + [C K R x *-i u - Q x «] U 1 + [C 2 R X *-2 U - C K Q x «.-\\ U 2 + 
+ [Cl R X 2 U - Cl Q x3 ] U K - 2 + [Cl R xu ~ C\ Q x2 ] U K ^+ 
'4D { +[ClR u 2- K 2 Q xu ]U 1 U K - 1 +[-C 2 +1 Q u ]U 2 U K - 1 + 

+ [R u - C\ Q x ] U K + [-Cl +X Q u ] U'U K + 
+ Remainder, 



where the term Remainder denotes the remaining terms in the expansion of R K . 

We note that the formula (1411 is valid for n = 3, comparing with ( 1401 . with the convention that 
the terms U K ~ 2 and U^~ x vanish (they coincide with U 1 and U 2 ), and replacing the coefficient 
—C 2 +1 Q u = — C\ Q u = —6 Q u of the monomial U 2 [/ K_1 by —3 Q u , as it appears in J40i . 
The proof goes by a straightforward computation, applying the recursive definition of this partial 
formula. 



3.3. Computation of R K in the general case. Following the exact same scheme as in the case 
n = 1 we give the general partial formula for R K . We start with the first three families of 
coefficients R^ , R ki k) and R^ k ^ k ^ . Let S q be the Rronecker symbol, equal to 1 if p = q and 
to if p ^ q. More generally, the generalized Kronecker symbols are defined by 5p\ , '"'p k k '■= 
§<n S q2 ■ ■ ■ S qk 

U Pl U P2 Pk 

By convention, the indices j, i\, i 2 , . . ■ , i\ run in the set {1, ... , to}, the indices k, k\, k%, ■ ■ ■ , 
k\ and I, l\, l 2 , ■ ■ ■ , l\ running in {1, ... , n}. Hence we will write Yl7^=i Y^T 2 =\ ' ' ' Y^,T=i as 

Y., v and Er i= i ElUi ••'ElUi as The letters i 1 ,i 2) ...,h and h,l 2 , . . . , l x 

will always be used for the summations in the development of R^ fc2 k ^ . We will always use 
the indices j and k% , k 2 , ■ ■ ■ , k\ to write the coefficient R^ k ^ k ^ . 
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We have: 



(42) 



ii Zi 

+ EE[-«e 



Zl z 2 



*i,*2 Zi,Z 2 



For R^, k2 we have: 



(43) 



E E [ 5 fc2 fi x fel «n + 5 fci fi x fe2 u*i _ 5 l &x kl , 



+ 



ii I 

Z^ Z^ "i2 



»li*2 Zi,Z 2 



)' 3 ■ 



Zl Z2 Z3 



+ E E 

11,12,43 h-h,h 

+ E E Kv.t ^4 - ^ ^1 - € <t q 1 x \] < 2 + 
+ E E Ktt^-4&^i-4tt^ 



il,l2 Zi,Z2, 



Since we also treat systems of order k > 3, it is necessary to compute fc2 fca . We write this 
as follows: 



(44) 



R£ 1 , fca , fc 3 = I + n + m - 



where the first term I involves only polynomials in : 



(45) < 



I = R> 

x ki x k 2 x k 3 



x k2 x ka u'i ^ u k 2 - Ll, x kl x k3 u t i k 3 x k± x k2 u'i 



I Ki x fc2 a; fc3 ti l i fc2 x 

ii d 

S j O h 1 TP 1 4- \ " U' 1 '' 2 ff- 7 + A' 1 '' 2 Z?- 7 + 



11,12 Zi,Z 2 

1 rd>(2 pi _ r'l i^ih _ rZi ^>Z 2 

i "%,fe 3 ""x^Ht^ °i 2 °fei ^x^x^un °i2°k 2 ^x kl 



x k3 u'i 



-5° S Ll O l2 

«2 &3 ^X kl X k2 V,'l 



U h U h + Z^ Z^ [%,fe2,fe3 u'lu^u'a' 

il,»2,»3 (1,(2,(3 

rj ^(1,(2 /oZ 3 _ ri rd,(2 W 3 

U i3 U k 1 ,k 2 ^x k ^u'iu t 2 U i3 U k 2 ,k 3 ^x, 



+ E E 



U,«2,«3j*4 ^1,^2,^3^4 



tl *2 *3 M ' 
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the second term II involves at least once the monomial W 1 , : 



(46) { 



h R j 



R 



2 Qx kl X k3 + Qx kl X k2 



+ / / / / L°fc 1 ,fc 2 ,fc 3 iiHu^ + °fcl,fe 2 ,fe 3 «n M -2 T °fcl,fc 2 ,fc 3 "llHl 



»ii*2 h,h,h 



, 5^2,23 Wl _r_ ci 2 ,Z3 Wl 

XfcjTt^ + °fc3,fcl ^X k2 U^ + %,fc 3 ^^11*2 



xi fx' 1 '' 2 n'3 n ( 3 + A' 1 '' 2 o' 3 + 



ti 12,13 



£ £ l ^3 (^fcl',fc2,fc3 < 3«*l„i 2 + ^fcl'.fa.fea < 5« i lt l «2 + 

jfj (<:h,l2,U Wi 



»li*2)*3 il,i2,i3.U 



u k u k 2 ,k 3 Via 



, ^(4,(3,(2 /q'i 



+ °fcl,fc 2 ,fc3 U il ^(2 U *3,*4 + Z^ Z^ [ °*2 ^°fcl,fc2,fe3 ^U 4 ! + 

21,22 ^1,^2,^3,U 

+a' 3 ' Zi '' 2 o' 4 4-a' 2 '' 3 ^ 1 "\1 Tp 1 TP 2 

^"kiMM ^ "k u k 2 ,k 3 ^ u *i ^i,i 2 ^3,(4 



and the third term III involves at least once the monomial Iff* l2 l;j (note that there is no term 
involving simultaneously l2 and l2 

111 ~ Z^ Z^ L°fel,fc2,fc3 «*1 °»1 V°fc2,fc3 ^fci + %,fcl ^fc 2 + 



(47) 



il il,( 2 ,(3 
+ ^ki,k 2 Qx k , 



tl,l 2 ,l3 



ci 2 ,i3j^4 yO/l 



«1,«2 (l,'2,i 3 ,i4 

A- 7 1 x' 1 '' 2 '' 3 n' 4 j- a' 4 '' 1 '' 2 n' 3 1 a' 3 '' 4 '' 1 n' 2 

°i2 \ °fci,fe 2 ,fe 3 ^u*l U k u k 2 M ^u'l ^ "kiMM ^u'l 



W'U! 2 , , . 
ti 12,13,14 



Before giving the partial expression of i? K we introduce some notations. For p G N with p > 1, 
let S p be the group of permutations of {1, 2, ... , p}. For q € N with 1 < g < |> — 1, let be 
the set of permutations a G Sj, such that cr(l) < er(2) < • • ■ < cr(g) and o-(q + 1) < a(q + 2) < 
• • ■ < <j{p). Its cardinal is C^. Let £ p be the group of cyclic permutations of {1,2, ... ,p}. 
Reasoning recursively from the formula of k2 k3 given by J44i . we may generalize Lemma|3] 

Lemma 4. For every n > 4 and for every j — 1, . . . , m, k\, . . . ,k K = 1, . . . , n, we have: 
(48) R fe 1 ,fe 2 ,...,fe re = 1" h + Remainder 



where h = R{ kiXk2 _ XkK , 



11 h 



ft 



r/'i 



o-eei 
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ii h,h 



< E C„ o 



(1) " :Cfc CT (2) "' Xk a( K ) 



U, 



h,h> 



*=E E 

ii h,...,l K -2 



E * 



ii, ,k 



R 



( 1 ),...,fc ff (»-2) ^(re-l^M*)"' 1 



< E 4 



Cl)l»M*^(»-8) ^ iE *V(,,_2) a! * 1 r(«-l) a!fc «rC<0 



ll,. 



*=£ E 

ii li,...,l K . 



E * 



d, 



R 3 



*=£ E 

21,22 ll,...,l K 



J,- d i ... .J 



Z^ %, °ii fc„ 



t££„ 



verged 1 



( )•' 

(l),...,fc„ C »_l) ^Xfc u<2 



Veer 1 V 

* = E E K (fc,t +i e 2 '+ : • + fa;' 2 <&) 



^ ^ 



Q i2 



»li»2 ^3, 



" i2 1 Z^ %, ,K ^un 

ree? 



/g = E E 



Tpl 



S h,...,ln n] _g] \ "* gh, W K 

E E K fct 1 - € (fct + • • • + fc.l ) 

11,12 ^lj"-)^K.+ l 



Li (2,...,t re +l 

anc/ where the term Remainder denotes the remaining terms in the expansion of R kl fc , k ■ 

In Iq the summation on the upper indices (li, . . . , l K ) gets on all the circular permutations of 
{1, 2, . . . , k} except the identity. In Ij the summation gets on all the circular permutations of 
{2,3,...,k + 1}. In Iq the summation gets on all the circular permutations of {1,2, ...,« + 
1} except the one transforming (1%, I2, ■ ■ . , l K +i) into (I2, 13, ■ ■ . , h). For k = 3, comparing 
with (1441 . we see that the formula remains valid, with the same conventions as in the case n = 1. 



3.4. Lie criterion and denning equations of ©tjm(£). We recall the Lie criterion, presented in 
Subsection 2.6 (see Theorem 2.71 of 1181 1: 
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A vector field X is an infinitesimal symmetry of the completely integrable system (£) if and 
only if its prolongation X^ of order k is tangent to the skeleton Ag in the jet space J^m- 

The set of infinitesimal symmetries of (£) forms a Lie algebra, since we have the relation 
[X.X 1 ]^ = [X^\X' (k) } (cfi ED). We will denote by 6tjm(£) this Lie algebra. The aim of 
the forecoming Section is to obtain precise bounds on the dimension of the Lie algebra &t]m(£ ) 
of infinitesimal symmetries of {£). For simplicity we start with the case n = m = 1. 

4. Optimal upper bound on dim K 6tjm(£) when n = m = l. 

4.1. Defining equations for 6rjm(£). Applying the Lie criterion, the tangency condition of 
to Ag is equivalent to the identity: 



(49) 



R' 



ox au 



R 



dF 

am 



OF 



R h 



dF 

dU K - 



= 0, 



on the subvariety Ag, that is to a formal identity in K{x,u, U , . . . , U K }, in which we 
replace the variable U K by F(x,u, U 1 , . . . , C/ K_1 ) in the two monomials U K and U 1 U K of 
R K , cf. Lemma [5] Expanding F and its partial derivatives in power series of the variables 
(U 1 , . . . , U K ~ l ) with analytic coefficients in (x, u), we may rewrite ( I49> as follows: 



(50) 



E 

..,/i K _i>0 



[^/il, ...,/*«_! u , (Qx k u')k+l<K, (R x k u')k+l<n)] x 

x(C/ 1 )^ 1 ...(C/ K - 1 )^- 1 =0, 



where the expressions 

(51) (ClVlfc + KK.t^'uOfc+Kli) 

are linear with respect to the partial derivatives ((Q K k u i)fc+i<rej (-^£c fc «0fe+(<«)' w ith analytic 
coefficients in (x,u). By construction these coefficients essentially depend on the expansion of 
F. The tangency condition d50t is equivalent to the following infinite linear system of partial 
differential equations, called defining equations of Srjm(£): 

(52) ...,/!„_! (x, u, {Q x k u i (x, it)Wz< K , {R x k u i (x, u)) k+ i< K ) = 0, 

satisfied by (Q(x, u), R(x, u)). The Lie method consists in studying the solutions of this linear 
system of partial differential equations. 

4.2. Homogeneous system. As mentioned in the introduction, we focus our attention on the 
case k > 3. Denote by (£q) the homogeneous equation = of order n. The general solution 
u — Yli=o ^ l • T ' cons i sts of polynomials of degree < k — 1 and the defining equation ( 149 \ 
reduces to R K = 0. Using the expression (I4H . expanding (15 0> . ( 15 1> and considering only the 
coefficients of the five monomials ct., U K ~ 2 , J7 re_1 , U 1 C/"" 1 and U 2 t/ K ~\ we obtain the five 
following partial differential equations, which are sufficient to determine ©rjm(£o) : 

R x * = 0, 

(«-2) 



(53) 



3 

(«-l) 



7?„ 



The general solution of this system is evidently: 

J Q = i + Bi + Ci 2 , 
[R=(k- l)Cxu + Du + E° 



(54) 



Q x 3 = 0, 
Q x 2 = 0, 

Q X u 0i 
Q„ = 0. 
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where the (n + 4) constants A, B, C, D, E°, E 1 , . . . , E^^ 1 are arbitrary. Computing ex- 
plicitely the flows of the (k + 4) generators d/dx, xd/dx, x 2 d/dx + (k — 1) xu d/du, ud/du, 
d/du, x d/du, . . . , d/du, we check easily that they stabilize the graphs of polynomials of 
degree < k — 1. Moreover they span a Lie algebra of dimension (k + 4) and the general form of 
a Lie symmetry is: 



4.3. Nonhomogeneous system. Consider for k > 3 the equation ( 1491 after replacing the vari- 
able U K by F. Let $(£/ A ) denote an arbitrary term of the form <j)(x,u)U x , where <j)(x,u) 
is an analytic function. We consider the five following terms <£>(ct.), <J>([/ K ~ 2 ), $(J7 K_1 ), 
^(U 1 and <J?(t/ 2 U K ^ 1 ). Since some multiplications of monomials appear in the ex- 

pression (@9}, we must be aware of the fact that <&{U 1 U K ~ 1 ) = ^(E/ 1 ) <I>([/ K-1 ) and 
$((7 2 t/ K_1 ) = <J>([/ 2 ) $([/ K ~ 1 ). Consequently in the expansion of d49l we must take into ac- 
count the seven types of monomials $(ct.), <£>([/ 2 ), <£>([/ K ~ 2 ), $(C/ K_1 ), ^(U 1 U^ 1 ) 
and $((7 2 [Z^ 1 ). The derivatives <9F/<9x, aF/9u, dF/dU 1 , dF/dU^ 1 appearing 
in the brackets of ( 1491 , and the term F appearing in the expression of R K after replacing U K by 
F (cf. the last two monomials U K and U 1 U K in (14 1\ ) may all contain the seven monomials ct., 
U l , U 2 , U K ~ 2 , U K ~ l , U 1 U K ~ l and U 2 U K ~ l . For F and its (k + 1) first derivatives we use the 
generic simplified notation 



(56) $(ct.) + ^(U 1 ) + <P(U 2 ) + $(J7 K ~ 2 ) + $(t/ K_1 ) + ^(L/ 1 C/ K_1 ) + <S>(U 2 U** 1 ), 



to name the seven monomials appearing a priori. Hence, expanding ( I49K picking up the only 
terms which may contain the five monomials we are interested in, and using the formula of 
Lemma[3]for R A (1 < A < k), we obtain the following expression: 



(57) 

' tlx* + [C 2 K R X 2 U - Cl Q x s] U K - 2 + [Cl R xu - C\ Q x i] U K ~ X + 
+ [C\ R U 2 - k 2 Q xu ] U 1 U*- 1 + [-C 2 K+1 Q u ] U 2 U K ~ l + 
+ {R u -ClQ x + [-Cl +l Q u ] U x }x 

x |$(ct.) + ^(U 1 ) + <S>{U 2 ) + <5>(U K - 2 ) + $({7 K_1 ) + ${U l U R ~ X ) + $(C/ 2 U^ 1 )} - 
- {Q + R + R x + [Ru - Q x ] U 1 + R X 2 + [2R XU - Q x =] U 1 + 
+ [R u - 2Q X ] U 2 + --- + R x «-3 + [Cl_ 3 R x ,-, u - Q x ,- 3 ] U x + 
+ [C K -3 R x *- 5 u ~ C re -3 Qx"- 4 -] U 2 + R x *-2 + 
' + [Cft-2 Rx^-^u ~ Qx"- 2 ] U 1 + [C«-2 R X "- i u ~ Qx"" 3 ] U 2 + 

+ [Ru - C\_ 2 Q x ) U K - 2 + [-<£_! Q u ] U 1 U K - 2 + 

+R x n-1 + [C K _ X R X *--2 U - Q x "-l] U + 

+ \Pk-\ Rx*-3 U - C re -i Qx"- 2 ] U 2 + [C K _x R xu - C 2 _ l Q x 2] U K ~ 2 + 
+ [C X K _ X R U 2 - (k - l) 2 Q xu ] U 1 U K - 2 + [R u - C\_ x Q x ] U K - X + 
+ [-ClQ u ] tfU"- 1 }* 

x |$(ct.) + ^(U 1 ) + $(t/ 2 ) + $([/ K ~ 2 ) + $(C/ K " 1 ) + ^(U 1 U^ 1 ) + $(C/ 2 U^ 1 )} 



Here the term Remainder consists of the monomials, in the jet variables, different from the five 
ones we are concerned with. The first four lines before the sign "— " develop R K and the third 
line consists of the factor F replaced by (I56> . In the last line (note that this is multiplied by the 



(55) 




+ Remainder = 0. 
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nine preceding lines) we replaced the (k + 1) first partial derivatives of F appearing in ( I49> by 
the term i56\ which we factorized. 

By expanding the product appearing in this expression J57> . and equaling to zero the co- 
efficients of the five monomials ct., U K ~ 2 , U K ~ 1 , U 1 U K ~ 1 and U 2 U K ~ X , we obtain the five 
following partial differential equations 



(58) 







= IL(x 


Q i Qx 1 R-j Rx 1 • • • 1 Rx K_1 1 J 


^ K I ^X 2 U 


- C\ Q x 3 


= IL(x 


^? Q i Qx i Qx 2 i R-i Rx ) • * * ; R x K -i ) Ru i Rxu) 


C K R X U 


- C 2 Q x 2 


= IL(x 


n i Q: Qx-i R-, Rx, ■ ■ ■ , i?a;«-l, Ru), 


C K R U 2 


^ Qxu 


= U(x 


Q -, Qx? • • • 3 Qx k—i i Qu-i R, Rx, • ■ • R x k ~^ •> 



Qu 



= Il(x 1 u,Q 1 Q x , 



Ru-, Rxu, ■ ■ • •> R x k ~ 2 

, Q x k-2 , i?, Rxi • ■ • Rx* — 1 3 

j RxUi ■ • * ) Rx K ~ 3 u) • 



(59) 



Here by convention II denotes any linear quantity in Q, R and some of their derivatives, of the 
form 

II(a;, u, Q x a lu b 1 , . . . , Q x a p u b p > Rx c iu d i i ■ ■ ■ , Rx c iu d i ) = 

where <f>i and are analytic in (x, u). For instance, the differentiation of II(a;, u, Q, R, R u ) 
with respect to x gives the expression Yl(x,u,Q,Q x , R, R x , Rxu)- Let us intro- 
duce the following collection of (k + 4) partial derivatives of (Q, R) defined by 
J := (Q,Q x ,Qx 2 , R, Rx, ■ ■ ■ , Rx*- 1 > Ru)- The aim is now to make linear substitutions 
on the system (I58> to obtain the system d68l where the five second members depend only on the 
collection J. The desired estimate dirriK 6rjm(£ ) < k + 4 will follow from d68i . 

Let us differentiate the third equation of (I58> with respect to x. Dividing by C\ we obtain: 

(k-1) 



(60) 



R, r 



■ Q x 3 = U(x, u, Q, Q x , Qx 2 , R, Rx 



. Rx n i Ru i Rx 



(61) 



Solving R X 2 U and Q x s by the second equality in ( I58> and by d4.3l we find 

{Qx 3 — n(x, w, Q 7 Qx-, Qx 2 ■) Rx-j • ■ • i Rx K i Ru-, Rxu)? 
Rx 2 u = n(x,u, q, Qx, Qx 2 , R, Rx, ■ ■ ■ , Rx^ , Ru, Rxu)- 

Replacing R x k by its value given by the first equality in (1581 we obtain for Q x 3 : 

(62) Q x 3 = U(x, U, Q, Q x , Qx 2 ,R, Rx, ■ ■ ■ , R x "~ 1 ,Ru, Rxu)- 

If we write the third equality in d58l as 

(63) R xu — n(a;, u, Q, Q X ,Q X 2 ,R, Rx, ■ ■ ■ , Rx"- 1 , Ru), 

we may replace R xu in ( 162k This gives the desired dependence of Q x 3 on the collection J: 

(64) Q x 3 = II(x, u, 

Q, Qx, Qx 2 , R, Rx, - - - , Rx"- 1 , Ru)- 

We may now differentiate the equalities ( 1631 and ( I64t with respect to x up to the order I. At 
each differentiation we replace Q x 3, R xu and R X " by their values in d63i . in (I64> and in the first 
equality in d58i respectively. We obtain for I G N: 



(65) 



Q x i + 3 — n(x, ?/, Q : Qxi Qx 2 1 ^3 Rx? ■ • ■ i Rx k ~ i 3 Ru), 

7/, Qxi Qx 2 i ^3 Rx-i ■ • ■ i Rx*- 1 3 Ru)- 



J ^x lJrl u 
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Replacing these values in the fifth equality of ( I58> . we obtain 

(66) Q u = U(x,u, 

Q: Qx, Qx 2 ) R? Rx, ■ ■ • 7 Rx K_1 ) Ru)- 

By replacing the fourth equality of J58i we obtain finally 

(67) R u 2 — n(x, u, Q, Q x , Q x 2, R, R x , . . . ,R x ^-i,R u ). 

To summarize, using the first equality of ( I58K using ( I66> . ( I67K d63l and ( I64l i. we obtained the 
desired system: 

R X K U, Q , Q x ) Qx 2 1 R^ F^X , • ■ • , Rx* — 1 ) -^U ) ; 

— n(x, li, Qa: ; Qx 2 1 R: Rx J • ■ • ) R X K—1 , -^U ) ; 

(68) < i?«2 =11(0;, u,Q,Q x , Q X 2, R,R x ,...,R x *-i,R u ), 

Rxu — U, Q , Qx , Qx 2 1 R> ^x , • ■ • j R x k—1 , Ru ) ; 

k — u, Q , , (5a;2 , R x , . . . , R x k—i , ) . 

We recall that the terms IT are linear expressions of the form d59t . Let us differentiate every 
equation of system (1681 with respect to a; at an arbitrary order and let us replace in the right 
hand side the terms R x *., R xu and Q x 3 that may appear at each step by their value in d68l . 
and then differentiate with respect to u at an arbitrary order. We deduce that all the partial 
derivatives of the five functions R x k, Q u , R u 2, R xu and Q x 3 are also linear functions of the 
(k + 4) partial derivatives (Q, Q x , Q x i , R,R X , . . . , R x *-i, R u ). Thus the analytic functions 
Q and R are determined uniquely by the value at the origin of the (k + 4) partial derivatives 
(Q, Qx, Q x 2 , R, Rx, ■ ■ ■ , Rx*- 1 , Ru)- This ends the proof of the inequality diiriK Srjm {£) < 
k + 4. □ 



5. Optimal upper bound on dim K 6t)tn(£ ) in the general dimensional case 



5.1. Defining equations for &t]m(£ ). In the general dimensional case, the tangency condition 
of the prolongation X K of X to the skeleton gives the following equations for j = 1, . . . , m and 

k\, . . . , k K = 1, . . . , n: 



(69) 



^-fci,...,fc„ 



.1 = 1 



dxi 



EE R * 



dFi 

n fci 



dUi 1 
'1 



<>'i /•-. 



E 



E 



IV, 1 



1 dU! 1 
ii,. 



on Ag, by replacing the variables l by ly* l wherever they appear. Let us expand 

Fjl i k and their partial derivatives and use the fact that R*^ k ^ are polynomials expressions 
of the jets variables (U{* , . . . , U£ ^), with coefficients being linear expressions of the partial 
derivatives of order < A+l of Q l and BP. We obtainfor j = 1, . . . , m and k\, . . . , k K = 1, . . . , n 
some identities of the form 
(70) 

Z2 ®ku'..,k K ; h, f^j") (Q l xOL u p)l<l< n ,\ a \ + \0\<K+l, (Ric, u 0)l<3<m, |a| + |,3| <K+1 ) X 



. U^ 1 X U; 



■u 



-1 Vl+2f.2 



0. 



satisfied if and only if the functions Q l and Ri are solutions of the following system of partial 
differential equations 



(7i) k:: 



(x, U, (Q x <x u 0)l<l< n , | Q | + |/3|< K +1, (^o u ,a )l<j<m : |a| + |/3|<K + l) - 0. 
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5.2. Homogeneous system. We start by giving the general form of the symmetries of the 
homogeneous system in the case n = 2. Then we prove the equality dimK(6r)m(£o)) = 



n 2 + 2n + m 2 + m C^_ l _l_ 1 in the case k > 3. 



In the case k 



n+K- 

2 we obtain: 



(72) 



l (x,u)=A l + B l kl x kl 

fci=l 



E 

ii=l 



+ V D kl xi x kl + V E h xi u n 



fci=i 



«i=i 



fci=l 



E 

ii=l 



my 



^ -Dfci a^fej + 
fei=i 



m 

E^ 

»i=i 



Here the (n+m)(n- 



-2) constants A 1 , B l 



feii 



C^, Dfe, F J , G J fei , fl? G IK are arbitrary. 



Moreover one can check that the vector space spanned by the (n + m) (n • 



m - 



2) vector fields 



(73) 



a 



dx 



fei 



s 



d 



x kl Xi 



dx k2 

d 
dx\ 
d 

Xi ~ h • 



+ x 



dx kl 

d 

dx n 

d 
dx n 
d 



+ u x 



d 
d 



+ u" 



d 
d 



x kl 



du* 1 



is stable under the Lie bracket action and that the flow of each of these generators is a Lie 
symmetry of the system (So). This proves that Srjm(£o) is indeed a Lie algebra with dimension 
(n + rn)(n + m + 2). Finally the corresponding transformations close to the identity mapping 
are projective, represented by the formula: 



(74) 



(x,u) 



k=1 «z,fc x k + 2_j i=1 ai,n+i u 



V 1 + ELl Tfe X k + Yh=1 In+i U* 

1 + ELi 7fc x k + YT=i Tn+i u z 



KKr 



l<j<m, 



0, the graphs 



It is clear that these transformations preserve all the solutions of (So) 
of affine maps from K™ to K m . 

In the case k > 3 we consider the homogeneous system (So) in which the second members 
Fl k vanish identically. Its solutions are the graphs of polynomial maps of degree < (k — 
1) from K™ to K m . The defining equations of its Lie algebra of infinitesimal symmetries are 
Rj. k =0, after having replaced the variables l by = l in ig and Iq in J48i . 
We will keep in this system the only equations coming from the vanishing of the coefficients of 



the five families of monomials ct., XJ\ 1 l 



2 '1, 



h h,-,l. 



andt// 1 , W 2 , (this 



is inspired from the computations in Subsection l4.2> . The coefficients of these five monomials 
families already appear in the expression ( 148 \ . Moreover we fix l\ = I2 = ■ ■ ■ = l K +i = I an d 
ii = 12, except for the coefficient of the monomial Iff 1 XJ\ 2 where we fix first i\ — 12 and 
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then i\ ^ %i. This provides the six partial differential linear equations: 
f = Ri , 



k-i ko k K ' 



°= E t:;;: 



1 i? 3 



-V E fc;::' 



— >«V( K -3) " a:fc CT ( K -2) a:fc CT (» C -l) ;r ' c <T(fc) 



L J *tl(l) fi <l(l(-l) ^^CTf^) 



(75) I 



-S j \ V ,l O l 

"ill u fc CT (i),-,fe CT (»-2) ^^.(..D^,,,, 

\<reer 2 

ii /,-< { r?u i, — ks{ [ y 4' 



Vo-eelT 1 



Veer 1 / 



kS 1 ( E fcT-ic--!)^ 



[0 = -c^ti^' 



To solve the system l !75> we fix the indices ki = • • • = k K = I and j = i\ in the sixth equation, 
implying Q l u , 1 = 0. Hence the terms following 5f i and Sf 2 in the fourth and in the fifth equations 
vanish identically. Let us choose the indices k\ = ■ ■ ■ = k K in the fourth and the fifth equations 
(this last equation is satisfied only for i\ ^ ii). We obtain first three simple equations, without 
any restriction on the indices: 



(76) 



= R? 



Finally we specify the indices in the third equation of i75i as follows: I = k K = ■ ■ ■ = k% = 
k-2 = ki ; then I = k K = ■ ■ ■ = &3 = ^ k\ ; finally I = k K = ■ ■ ■ = k%, k% ^ k2, k% ^ k\. 
This gives the three following equations: 



(77) 



We specify the indices in the second equation of H51 as follows: / = fc K = • ■ ■ = £3 = k^ = k\ ; 
then I = k K = ■ ■ ■ = £3 = k^ 7^ k\ ; then / = fe re = ■ • • = ks, k% ^ k^, k% 7^ fci ; finally 
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(78) 



= ki, I 7^ ki, I 7^ k%, I 7^ k%. This gives the four following equalities: 
( = C 2 K R 3 

= C*_! i? J 4l - C*_j # Q* 2 t Xh Xh , k 2 + ki 



C 3 ft 3 O kl 







I ^ h, Ij^ki, l^ks 



kaj^ki, k 3 ^ k 2 , 



Let us differentiate now the equations J77> with respect to the variables xi as follows: we dif- 
ferentiate J77> i with respect to ie^ ; then we differentiate J77b with respect to £fc 2 ; finally we 
differentiate J77h with respect to Xk 3 - This gives the three following equations: 



CI 2 i) j (l kl 



(79) 



= Cl R 3 

x kl x kl u l i 

= R 3 h - Cl , Q^ 2 _ _ , fc 2 ^ fci, 



(5? O ka 

»i ^ fcl a; fc2 :E fc3 ' 



&3 7^ fel, ^3 7^ fc 2- 



The seven equations given by the systems (I78> and ( I79> may be considered as three systems of 
two equations (of two variables) with a nonzero determinant, to which we add the last equation 
J78t zi. We get immediately: 



(80) 



'0 


= R j 

x kl x kl u l i 


- S 3 O kl 

n ^x kl x kl x kl i 






< 


= R j 

x k± x k2 u l i 


— Qx 2 kl x k2 x k2 7 ^2 







= R j 

x kl x k2 u z i 


- S 3 O ka ki 

u il ^x kl x k2 x k3 -i "-3 


7^ ki , fc 3 7^ fc 2 


^ 


- 6 j O l 

11 ^x kl x k2 




l + k z . 



It follows from these relations and from the relations Q l u — R J u io = obtained in M6\ that 
all the third order partial derivatives of Q l vanish identically, this being also satisfied by the third 
order partial derivatives of R 3 containing at least one partial derivative with respect to u 11 : 



(81) 



= O l = O l = O l ■ ■ — o l 

U ^X kl X k2 X k3 ^X kl X k2 U z l ^X kl U'lu'2 ^Ii'lB'2u'3l 



= R 3 

X k , I t „ u'l 



R 3 



kl U i lU i 2 



R 3 

U*lU»2u»3 



It follows from the equations J76i and dH 1 b that all the functions Q l are polynomials of degree 

< 2 with respect to the variables Xk ± and all the functions R 3 are a sum of a polynomial of degree 

< (k — 1) in the variables Xk ± and of monomials of the form u 11 and Xk x u 11 . Let us develop now 
the relations 1771 separately for j = i\ and j ^ i\. We obtain the five equations: 



(82) 



(0 = 01 R' 1 u - c 2 K Q k x \ Xk , 

Q =K\ lU ^- c UQw k2 ^ k ^ 







k 1: k 3 7^ fc 2 
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According to the equations ( I76> . ( I8U . ( I82> . we have the following form of the general solution: 

n n 

3 z (x,m) = A 1 + ^2 B k x x k! + Gk i Xfc i Xh 



(83) I 



fci=i 



fci=i 



fci=i 



= ^ (k - 1) C kx x kl u 3 + K uti + Ej '° + Yl E H x ^ 

k\ — 1 i± — 1 

+•■■+ E 

l<fci<-<fc K _i<n 



£ 'fel,...,fe«_l ' ' ' ^fc^-l ' 



Here the n + n 2 + n + m 2 + mC^ l+ l_ 1 constants A 1 , B l ki , C kl , D J h , E j <°, E^f, 
Ej.'*~ k 6 K are arbitrary. Moreover one can check that the vector space spanned by the 



vector fields 



(84) 



d 



d 



dx kl dx k2 

a 



Xk > [ Xl + 

8 d 



■>n— + («-!) U 



OX n 




du r - 



^ du 12 ' du 11 ' kl du n 



x kl ---x kK _ ± dun , 



is stable under the Lie bracket action and that the flow of each of these generators is indeed a Lie 
symmetry of the system (So). Finally the Lie symmetries of (So) have the following form: 
(85) 

) )■ 

' l<j<m J 

0, namely 



{ [1 + E*=i SkXk]*- 1 

We note again that these transformations preserve the solutions of (So) : u% k ,.. Xk = 
the graphs of polynomial maps of degree < (k — 1) from K" to K m . 



5.3. Nonhomogeneous system. Let k > 3. Let us expand the defining equations ( I69i as 
done in ( I70> . We will write only the coefficients of the five monomial families ct., l , 

Ufl t^, U\\ U t,...,i K and u il.i 2 U it i K j- Moreover, we fix always Z x = h = • • ■ = /« = 
/ K+ i = I and ii = 12, except for the fourth family of monomials where we distinguish the 
two cases i\ = %2 and i\ ^ 1%. Thus we obtain six linear equations of partial derivatives, the 
members on the left side (coming from the expression of k given by Lemma |4} coin- 

cide with the members on the right hand side of (175 \ . Furthermore, the members on the right 
hand side are exactly the same as those obtained in J58i . with more indices! We use the letters 
I', k[, . . . , k' K — 1, . . . , n and j', i[ = 1, . . . , m for the indices of the arguments of the expres- 
sions II, obtaining the six following equations, which generalize the equations J58l i: 



(86) [1] 



= H lx,u,Q l ,Qi ,JP ,R*. 



•Ml' 
U 1 
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M : R3 x kl x k2 -x kK - 11 \ Q L , Qx k ^ RJ > R3 x kli ' ' ' ■ ,R ^ k ' 1 -^ k ' R _ 1 ,R3 u i'i 

[2] : V Si' \ R j 



(i)v.*<r(»-a) :rfe „( K -i) :rfe „( re )" n 



"ill \(1) ^(.-S)^ 1 *,!,-!) 1 ^,-!) 1 *^!,) 

\aesr 3 

n[x,u,Q l \Qi k ,,Qi k , Xk ,,W' ,Ri k ,,...,Ri k ,... Xkl ,R j ',,R j ' 4 , 



[3] : V ^ 



(l),.-->«<7(«-l) at 



u h \ Z-^ "^(l),-,fc CT (.-2) ' 1, »i„ ( ,_ 1) il„ ( , 1 



\<ree 



K — 2 



= n [x,u,Q l ,Q l x ,,R 3 ,R J X R 3 X ,... x , ,R\, 



Wer 1 

= II ^X, U, Q , Qx,./ I • ■ • 1 Qx k ^---X k i ' *3 u »'i ' ^ ' ' ' ' ' ' —X fc / 

, It , It ./,..., it ./ I . 

M 1 I t l U 1 X,./ U 1 / 



[5]: 2K^;;;-;^tf, lii , 2 I ^C ( i),'..'fc<,( K -i) ^* a(K) «' 2 

Wer 1 

Veer 1 / 



II ( X, U, Q , Q Xk , , ■ • ■ , Qx fe / ■••x t ,/ i > Qyi'i ' ^ ' R x k ,^ ) • • ■ > ^x fc / ■••x fc / 



"1 "ic-l " "1 "1 re-1 



U 1 X k f U 1 X k f ■■■X k / U 1 



[6]: -^ +1 ^fclCi = 

= II I X, U, Q ! Qx,./ I • • • 5 Qx,./ -"X^ ' ^ ' I • • • ) R X k l ---XU 



,R ., ,R R ., 

U 1 X t / It 1 ! Itl -I t l U 1 

K l K l 



we get the following Lemma: 
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Lemma 5. Let J denote the collection of n + n 2 + n + m , re _ 1 + m 2 partial derivatives 

(87) J : , 0*^ , <,,... , <, ..^ .//:•). 

After linear combinations on the system (1861 we obtain the following equations: 



■U(x 


u, J) 


= R j 




U(x 


u, J) 


— Q u { i > 




U(x 


u, J) 


= R j ■ 




< n(x 


u, J) 


= Q l 




n(x 


u, J) 


= # j 4l , 

x kl u'i ' 




U(x 


u, J) 


— O fel 


fci ^ k 2 , 




u, J) 




I ^ ki, I ^ k 2 



Moreover all the partial derivatives ( with respect to xi and u L ) up to order three of the coefficients 
Q l and i?- 7 of the vector field X G 6rjm(£) are of the form H(x, u, J). Hence every function Q l 
and W is uniquely determined by the values at the origin of the n + n 2 + n + m C^7 K _ 1 + m 2 
partial derivatives \87\ . This implies that dimK Srjm(£) < n 2 + 2n + m 2 + m G^7 re _ 1 . 

Proof. Since the second part of Lemma[5]is immediate let us establish only the identities (1881 . We 
first specify the indices in the equation (I86> p] as follows: I = k K = ■ ■ ■ = ks = k% = ki ; then 
I = fc K = • • • = fc 3 = k 2 ^ k\ ; and finally I = k K = ■ ■ ■ = fc 3 , fc 3 ^ k 2 , fc 3 ^ k\. This gives 
three equations whose members on the right hand side are the same as those in the equation (I77i 
and whose members on the left hand side are the same as those in the equation (I86> [^] : 
(89) 

C « R l kl u'i ~ C * S n Qxl^H < 
R 3 u -C£_i& Q k J x . , k 2 =£ki, 

x k 1 ul 1 fc l k 2 ' 

We remark that these three equations (after specialization of j = i\ or of j ^ i\ and after 
some easy linear combinations) provide directly the fifth, sixth and seventh equations of J88I . In 
particular we may replace the values of the partial derivatives R J ., and Q x . x with k[ ^ k' 2 

or I' 7^ k[, I' k' 2 appearing in the expressions IT of the second member of {86} [i] by their 
values just obtained from the fifth, the sixth and the seventh equations of d88t . This gives the first 
equation of (1881 . 

Then we specify the indices in J86> p] as follows: I = k K = ■ ■ ■ = fc 3 = k 2 = k\ ; then 
I = k K = ■ • • = &3 = &2 7^ kx ; then I = k K = ■ ■ ■ = fc 3 , fc 3 ^ k 2 , fc 3 ^ ki ; and finally 
I = k K = ■ ■ ■ = fci, I 7^ k\, I 7^ ki, I k3. This gives four equations, whose members on the 



1 ^X, I ■ 


T>j' 
K—1 U 1 




'"V ' R 
re— 1 u 1 


' RJx k[ ■ 


■x k , ,R\> 

K—1 U 1 
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right hand side are the same as those in ( I78> and the members on the left hand side are the same 
as those in (j86j[ 2 ]: 



(90) 



n x,u,q 



pj p3 p3 p3 p3 

3 X./ x.i 3 ^ 3 ^x.i j ■ • ■ 3 J^x.i ■■•x,/ 3 *l 3 ^ 
re l K l ^2 1 1 k — 1 u 1 x,,u 1 



!' r\l' rJ-' r>j' pi' 



J ! | / . tl . Q , Qx k / J Qx k / x k i 1 ^ i ~^ x k' ' ' ' ' ' 



X , f ■ ■ -X ,/ 1 * *■ ^/ 3 ^/ 

1 K — l U 1 x.i U 1 



fc 2 / fcl, 



II | a;,u, Q ! ,QL 



/ > ™6 x ./ x , / j 3 "a; , / 3 ■ • ■ 3 "a; j —33./ 3 i ; 3 
1 1 2 1 1 k — 1 it 1 a:./ w 1 

"l 

R 3 - Cl_zSl Q k J._ , ^3 / ki, k 3 / fc 2 , 



pj' pj' 



II | U, Q iQxt/'Qx./X./'R 3 ^S, / 3 ■ • ■ 3 / ■■•X, / > R > R ^ 

1 2 1 1 k — 1 u 1 x.i u 1 



l=£ki, I k 2 , I ^ k 3 - 



il ^ x k 1 Xk 2 x k 3 > 

Using the fifth, the sixth and the seventh equations of J88i just obtained, we may replace the 
partial derivatives R J ., and Q l with k[ ^ k 2 or V / k[, I 1 ^ k' 2 appearing in the 

expressions II of ( I90> . providing four new equations in which the arguments of II are the desired 
ones: (x, u, J), where J is defined in (I87> : 

f n(s, u, J) = C7^ i? J , - C\ 61 Q\\ xi xi , 
IL(x, u, J) = CU K kiXk2uil - CU 51 Q k x l iXk2Xk2 , k 2 + kx, 

[n( !C) u,J)=-<5f 1 Qi Jbi , 1,4*1, l±k % , l^k 3 . 



(91) 



Let us differentiate now the equations d89l with respect to the variables xi as follows: first 
we differentiate d89t i with respect to ; then we differentiate (I89i -? with respect to Xk 2 ; 
finally we differentiate ( I89h with respect to Xk 3 - The arguments in the expressions II in the 
equation d89t contain now the terms R 3 Xk , ... Xk , ; we replace them by their value given in the first 

equation of d88l already obtained. The arguments also contain the terms i? J , and Q l x 

with k[ ^ k' 2 or /' ^ k{, I' ^ k' 2 . We replace them by their value given by the fifth, the sixth 
and the seventh equations of (I88i . We obtain three new equations in which the arguments of the 
expressions IT are the desired ones: (a;, u, J), where J is defined in d87l : 



(92) 



n (x, u, J) = C\ Ri tl - Cl Si Q k x \ Xk Xh , 

n (a:, u, J) = R j M - Cl_ t d? Q k x \ _„ , k 2 ^ k u 

n (as, u, J) = — ^ Q x l iXk2 x k;j 3 fca^fci, fc 3 ^fc 2 . 



The seven equations d9 1 1 and d92i may be considered as three systems of two linear equations of 
two variables with a nonzero determinant, the seventh equation being the last equation in ( I91l l. 
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We immediately obtain: 



(93) 



'n(x 


U,J) 


= R J 

x kl x kl u'i 


- S j O kl 

u l! ^Cx kl X kl X kl 




U(x 


u, J) 


= R j 


- S j O k2 

ll ^oc kl x k2 x k2 


k 2 ^ h, 


" U(x 


u, J) 


= R j 

x kl x k2 u z i 


- s j o ks 

Jl ^x kl x k2 x k3 




U(x 


u, J) 


— Hi Qx kl x k2 


x k3 , h^h, 





giving the fourth equation in i 

It remains now to obtain the second and the third equations in (1881 . Let us write firstly 
equation J86> [k] with the choice of the indices j = ix, I = k% = ■ ■ ■ = k K . This gives the 
equation: 



Q 



(94) 



u'l n ^X, U, Q , Qx h ,^ i ■ ■ • j Qx h /^ ■■■x k t ^ i I& ' ) ^c fc / ! ■ ■ ■ ' ^Xj,/ ■■■x k , ^: 



W , . h" /»'•' 

u l i x,j u'l x,j---x,j u 1 



We observe first that the differentiation with respect to the variables xi of one of the expressions 
n(:r, u, J) remains an expression Il(x, u, J). Indeed we see from d87i that there appears, in the 
partial derivative J Xl , derivatives Q l x , x , with k[ ^ k' 2 or I' ^ k{, I' ^ k' 2 . We may replace 

k l k 1 

them by their value obtained in the sixth and the seventh equations of J88b . It also appears some 
derivatives Q l x , X/Xl (we replace them by their value obtained in the fourth equation of J88I T). 



1 2 3 

some derivatives R Xy ■ ■ -x k , (we replace them by their value obtained in the first equation of (I88» 
and some derivatives R 3 ., (we replace them by their value obtained in the fifth equation 

xu u z i 



of ( I88t ). Consequently we may write: 

(95) [U{x,u,J)] xl =n(x,u,J), 

It follows that any derivative with respect to xi (to any order) of the fourth and the fifth equations 
of J88t > provides expressions of the form Il(x, u, J). In other words for any integer A > 3 and 
any integer /i > 1 we have 



(96) 



U(x, U, J) — Q Xki x k2 x k3 ... Xki , 



IL(x, u, J)=R J 

v ' ' ' x k , ---Xk u'l 



We may replace then these values in the equation J94-I . replacing also the derivatives Q l x , x , 

k l k 2 

with k' x ^ k' 2 or I' ^ k[, I' ^ k' 2 by their values obtained in the sixth and the seventh equations 
of d 8 8 1 > - This gives the second equation of d88t >- 

We also remark that by a differentiation with respect to the variables xi, the second equation 
Q i t — n(x, u, J) just obtained implies, using d95l : 

(97) Il(x,u,J) = Q l XHuil . 

It remains finally to write d86t ^] first with the choice of indices I = ki = ■ ■ ■ = k K , j = i\ 
then with the choice of indices I = k\ = ■■ ■ = k R , j ^ i\. We also write d86t [^] first with the 
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choice of indices I = k\ = ■ ■ ■ = k K , j = i-2 then with the choice of indices I = k\ 
3 h, j 7^ i%- We obtain four new equations: 



(98) 



R\ tl -kQ 1 u =n [x 



J i 5 ^ ) Mfa:,/ 7 ■ ■ ■ t ^ x,/ ■•■x. 



it 1 X, / 11 1 X./---X 



T>3 T>3 f>3 
) ; -^Xj^f 5 ■ ■ ■ 5 

R j ' 



^,11 -,41 n[^,W,Q , Qa; , , . . ■ , Qx,/ ■■■X,i J Q i' 1 R 1 ^X, / ! ' ■ ' 5 ^ 



, It ; -ft XL ^/ 

It 1 li 1 3? fc / ■■■X^f U 1 



2iT 



kQJ 1 „ij = n ( q ! ,q' , . . . . 



U 1 fci «i 



pj pJ pJ 

it 1 n, / u 1 X , / ■ ■ -x ,/ u 1 

1 1 K — 2 



ii 7^ *2, 



3 



Using the equations of J8 8t > we already obtained (namely all except the second equation), us- 
ing {%J and @7j, we may simplify these four equations: 



(99) 



U(x,u,J) = R\ M 

n(x,u,j) = r\, „ 

n(x, u, j) = r\ 4 , 
n(x,u,j) = r 3 i, 



This gives the second equation of J8 8t > . completing the proof of Lemma[5]and consequently the 
proof of Theorem^ D 
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